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ABSTRACT 

The collapse of a non-collisional dark matter and the formation of pancake structures 
in the universe are investigated approximately. Collapse is described by a system of 
ordinary differential equations, in the model of a uniformly rotating, three-axis, uni¬ 
form density ellipsoid. Violent relaxation, mass, and angular momentum losses are 
taken into account phenomenologically. The formation of the equilibrium configura¬ 
tion, secular instability and the transition from a spheroid to a three-axis ellipsoid are 
investigated numerically and analytically in this dynamical model. 

Key words: dark matter - large-scale structure of Universe. 


1 INTRODUCTION 


According to modern cosmological ideas, most of the matter in the Universe is in the form of so-called cold dark matter, 
consisting of non-relativistic particles. The study of the formation of dark matter objects in the Universe is based on N-body 
simulations, which are very time consuming. In this situation, a simplified approach may become useful, because it allows one 
to investigate many different variants and to obtain some new principally important features of the problem, which could be 
lost or not visible during long numerical work. 

The systematic analysis of ellipsoidal figures of equilibrium is presented in the classical book of lcha,ndra,sekharl ( 1 
Different types of incompressible ellipsoids are investigated there on the basis of virial equations, as well as their dynamical 
and secular stability. In particular, the points of onset of the bar-mode dynamica l and s ecular instabilities of the Maclaurin 
spheroid, for its transition to the Jacobi ellipsoid, are found. iLvnden-Beli JlQBdl applied the Chandrasechar’s virial 

tensor method for a rotating, self-gravitating spheroid of pressure-free gas and showed the growth of non-axisymmetric 
perturbations during collapse. It was also discussed that the slowly shrinking Maclaurin spheroid will enter the Jacobi series if 
it shrinks slowly enough for the dissipati ve mechanisms to b e ope rative. There have been nu merical investigations of collapsing 
pressureless spheroids in the papers of |Lwidei^Belj|_^^6^^_md lUin. Mestel fc Shul 1^6^. The dynamics of a non-rotating 
sphere in two dimensions was considered in lLvnder^^li Jl97ijl , and it was shown that the pressure prevents the development 
of large-scale shape instabili ty if initially the gra vity is more than three-fifths pressure resisted. There is a wide-ranging review 
of this topic in the paper of lf.vnden-Relll ( 1 99(ji in memory of Chandrasekhar. 

Sub sequent investigations of rapidly rotating figures are connected with stars, and with large-scale structure of the Uni¬ 
verse. IniRosensteel fc HuvUl99lll the virial equations for rotating Riemann ellipsoids of incompressible fluid are demonstrated 
to form a Hamiltonian dynami c al syst em. There is a detailed description of the ellipsoid model of rotating stars in the papers 
of lUai. Rasio fc Shaoirol (Il99,ill994i]l^ . Using a variational principle, they derived and investigated the equations for the 
evolution of a compressible Riemann-S ellipsoid, incorporating viscous dissipation and gravitational radiation. The solutions 
of these approximate equati ons permitted them to obtain equilibrium models, and to investigate their dynamical and secular 
stability. In a recent paper dsha.nirnlknnjh secular bar-mode instability driven by viscous dissipation was investigated, using 
these equations, and the point of instability of compressible Maclaurin spher oids was found. 

The modern theory of a large-scale structure is based on the ideas of IZerdovid] (Il97dl . concerning the formation of 
strongly non-spherical structures during non-linear stages of the development of gravitational instability, known as ‘Zel’dovich’s 
pancakes’. Numerical simulations for these objects have been performed subsequently by many groups. In the case of structures 
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in dark matter, we are dealing with non-collisional non-relativistic particles, interacting only by gravitation. The development 
of gravitational instability and collapse in dark matter do not lead to any shock formation or radi ation losses, but ar e 
characterized by non-collisional relaxation. This relaxation is based on the idea of a ‘violent relaxation’ of iLvnden-BelJ (Il967ll . 

Here we derive and solve the equations for the dynamical behaviour in a simple model of a compressible uniformly rotating 
ellipsoid. We derive the equations for axes with the help of variation of the Lagrange function of an ellipsoid. Correct description 
of pressure effects, attained by such an approach, and the addition of relaxation permit us to obtain the dynamics of motion 
without any numerical singularities. In our model, motion along three axes takes place in the gravitational field of a uniform 
density ellipsoid, with account of the isotropic pressure, represented in an approximate non-differential way. Relaxation leads 
to a transformation of the kinetic energy of ordered motion into kinetic energy of chaotic motion, and to increases in the 
effective pressure and thermal energy. All losses are connected with runaway particles. The collapse of the rotating three- 
axis ellipsoid is approximated by a system of ordinary differential equations, where relaxation and losses of energy, mass 
and angular momentum are taken into account phenomenologically. The system is solv ed numerically for s evera l parameters, 
characterizing the configuration. The approach in this work is similar to that used in iBisnovatvi-KogarJ 1 I 2 OO 4 I 1 . where only 
dark matter spheroids (a = 6 yf c) were considered. In this case there are analytical formulae for the grav itational potential 
and fo rces. In the description of violent relaxation and different kinds of losses, we use the same approach as iBisnova.tvi-KoganI 
ll2nn4l . In the present paper we also find the point of the onset of secular instability of a compressible Maclaurin spheroid, 
using the derived dynamical equations, and also from the analysis of the sequence of equilibrium configurations. A simple 
analytical formula for the point of onset of the bar-mode secular instability of the Maclaurin spheroid is found. 


2 EQUATIONS OF MOTION 

Let us consider a three-axis ellipsoid, consisting of non-collisional, non-relativistic particles, with semi-axes a yf 6 yf c 
z'^ 

— + TJ H- 2 ~ 

and rotating uniformly with angular velocity fl around the a-axis. Let us approximate the density of the matter, p, in the 
ellipsoid as uniform. 

The mass m and total angular momentum M of the uniform ellipsoid are connected with uniform density, angular velocity 
and semi-axes as {V is the volume of the ellipsoid) 


= pU = ^ pabc, M = pfl {x^ + y^) dV = ^ + b^) . 


( 2 ) 


We assume a linear dependence of the velocities on the coordinates in the rotating frame. 


by 


cz 

c 


The gravitational energy of the uniform ellipsoid is defined as llLandau fc Lifshit jllQQ.lfl 

00 

3Grrd f du 


( 3 ) 


( 4 ) 


10 J G (a^ -I- u){W' + u)(c^ + u) 

0 

and is expressed in elliptical integrals llOradshtein fc Rvzhiklir962h . Let us consider a compressible ellipsoid with a constant 
mass and angular momentum, a total thermal energy of non-relativistic dark matter particles Eth ~ ~ and 

the relation between pressure P and thermal energy Eth as Eth = In absence of any dissipation, this ellipsoid is a 

conservative system. To derive the equations of motion, let us write the Lagrange function of the ellipsoid: 


L — Ukin Upot , Upot — Ug -f Eth Urot , 


Ukin = P {Vx + Vy + v^) dV = — {d -1-6 -be) 


Eth = 


Eth.ini^^inbinCin) ^ 

(abc)^/^ 




where the entropy function e is constant in the conservative case, but variable in the presence of dissipation, and 
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Urot ^ J ^rotdV ^ ^pQ^ J {x^ +y^) dxdydz = ^ -b Y) . 

V V' 

Taking into account the second expression in 0, we obtain the relation 
5 
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By variation of the Lagrange function we obtain Lagrange equations of motion: 

OO 

du 


“iGm 

a = --— a 


10 1 £ 25M2 a 

+ ^ - . , xo/, + ■ 


2 J (o^ + u)A 3m a {abc)'^/^ nfi (a^ + ’ 

0 


(9) 


b = - 


3Gm 


du 


10 1 £ 


25M^ 


{b'^+u)A 3m b (abc)^/^ m? (a^ + 6^)2 


( 10 ) 


3Gm 

c =--— c 


du 


+ 


10 1 £ 


(c^+u)A 3m c (abc)'^/^ 


A^ = {a + u)(J) + u){c + u) . 


( 11 ) 


It is easy to chec k that the equilibrium solution of these equations are the Maclaurin spheroid and the Jacobi ellipsoid 
llChandrasekhai^ll969^ . 


3 EQUATIONS OF MOTION WITH DISSIPATION 


In reality there is relaxation in the collisionless system, connected with phase mixing, which is called ‘violent relaxation’ 
llLvnden-BelJll96^ . This relaxation leads to a transformation of the kinetic energy of ordered motion into kinetic energy of 
chaotic motion and increases the effective pressure and thermal energy. 

The main transport process is an effective bulk viscosity. Therefore, there is a drag force, which is described phenomeno¬ 
logically by adding the terms 


a b c 

'l~ret '^rel '^rel 

on the right-hand sides of the equations of motions. 

Here we have scaled the relaxation time Trei by the Jeans characteristic time 


T - - 


27r 


= 27r 


abc 


V AnGm V 3Gm 

with a constant value of arei, 


'^rel — 27r Otrel 


abc 

3Gm 


( 12 ) 


(13) 


(14) 


The process of relaxation is also accompanied by total energy Utot, mass and angular momentum losses from the system. 
These losses may be described phenomenologically by characteristic times Tei, Tmi and tmi, scaled by rj with the coefficients 
del, dmi and omi: 

Utot ~ Ukin T Gg J- Urot T d^th- 

The corresponding equa tions describing the various losses are derived for a three-axis ellipsoid in section 4 of the paper by 
iBisnovatvi-KoganI J2004ll . 

Using the variational principle and derivi ng the equation for the entropy function from energy balance give the correct 
relations for the pressure and the total energy. iFuiimotol il968h derived the equations of motion for a rotating ellipsoid from 
a system of hydrodynamical equations. However, the accounting for pressure effects in his approach was inconsistent with 
thermal process es, l eading to the wrong resu lts for the dynamical behaviour of a system with radiative losses (compare with 
iLai et aHll994al and iBisnovatvi-KogarJDOOdh . 


4 NON-DIMENSIONAL EQUATIONS 

To obtain a numerical solution of the equations, we write them in non-dimensional variables. Let us introduce the variables 

P 


a ~b,c , m ~ M , p ~ U - Eth 
T ’ b = , c= — , m = — , M = — , p = — , U = — , Eth = -r- 

to ao ao ao mo Mo po uo Uo 

The scaling parameters to, ao, mo. Mo, po, Uo, Ho and £o are connected by the following relations: 

Mo 


£0 


IL 

to 


3 

,2 _ “O 
to — 


Gmo 


Mo = Gaoml, 


Uo = 


ao 


mo 

Po ~ TT’ “o - 


- 9 , £0 — Uoao- 
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The parameter Uo is used for scaling of all kinds of energies. Hereafter the ‘tilde’ is omitted for simplicity. In non-dimensional 
variables, we have 

47r TTl „ . 2 t2\ 

m=—pabc, M=—D,{a+b), Ti = 2n ai 
o 0 

Taking into account violent relaxation, total energy, mass and angular momentum losses, the dynamics of the system is 
described by the following non-dimensional system of equations: 



a 


a 


.. _ a dm 3m j du 

m dt 2 J (o^ -I- it)A ' 3m a (abc)'^/^ ' m^ (a^ -|- b'^)'^ ’’'rei ’ 

0 


10 1 e 25M^ 

+ - - , . + 


(15) 


- b dm 3m 


du 


10 1 e 25M^ 

+ 7^ T . , + 


m dt 2 7 ( 6 ^ + u)A 3m b (abc)^/^ m^ (a^ -|- 52)2 ’ 

0 


(16) 


.. _ c dm 3m 


du 


+ 


10 1 


m dt 27 (c^-I-u)A 3m c (abc)^/^ ’ 

0 


£ = (abc)^^^ Uk 


\( ^ 

1 

2 ) 

Urot 

( 2 

1 ^ 

Ukin 

UgTml_ 

V T-pei 

Tel 

'^ml J 

U, 

\tmi 

'^ml J 


(17) 

(18) 


1 . 2 , ^2 , . 2 ^ / / I du \ M 

+■=*4 / A I' M- 


{—)' 

\minJ 


, = {a + u){l} + u){c + u ). 


(19) 


This system is solved numerically for several initial parameters. 


5 NUMERICAL RESULTS 

The integrals describing gravitational potentials and forces, are expressed in elliptical integrals llOradshtein fc Rvzhikll^962^ . 
For ellipsoids that are close to a spheroid or to a sphere, we have carried out an expansion to linear terms near points a = b ^ c, 
b = c^a, c = a^b and a = b = c, respectively. These expansions include only analytical expressions. All integrals reduce to 
two, 10 and Jl, i.e. 


70 = 


du 




71 = 


du 


{a? + u) \J {a? + u) (K + w) (c^ -I- u) 


0 0 

In the limiting cases these integrals are expressed analytically as follows: 

arccos (c/a) 


,,_,arccos(c/a) 1^,2 „ 2 ^ 
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71 = - 
at 

^-1 

a 


( - 
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{c/a) ^ 

\ 1 

_ 1 
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1 c 

3 

arccos(c/a) 

\^a2 


-c2 , 

/ 

2 

4 a2 (a2 — c2)^ 

2 0“ (a2 - 

c2) 4 

(a2_c2)5/2_ 
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71 = - 
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a 
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a 
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Table 1. Initial values and parameters of the calculated variants. 


Cases 

1 

2 

3 

4 

Figure 

1,2,3,4 

5 

6 

7 

a 

1 

1 

1 

1 

b 

1 

1 

1 -1- 10-® 

1 -1- lo-'i 

c 

1 

1 

1 + 2 X 10-® 

0.5 

^in 

0.2 

0.5 

0 

0 

^in 

0.2-1-10-1° 

0.5-1-10-1° 

0 

0 

^in 

0.2 

0.5 

0 

0 

^in 

0.5 

0.1 

0 

0.317 

^in 

0.01 

0 

0.01 

0.149 

^rel 

3 

3 

- 

- 


15 

15 

- 

- 


at 


- 1 


71 = 2 


at 

^-1 

c 

70= - 
a 

at 


< 0 . 0001 , 


cosh 


b,c> a 
{a/b) 


j b"^ — 


- 1 


^ 0 . 0001 : 


(23) 




■ + 2 


- 3 


cosh ^ (a/b) 


(b^-a^)^ a (62- a 2)2 (^2 _ ^ 2 ) 5/2 


< 0.0001, b,c < a, 


^-1 

a 


^ 0 . 0001 : 


(24) 


1 

+ 3 


a^-b^ 


71 = 


2 


+ 5' 


1 

+ 5 




- 1 


< 0 . 0001 , 


-1 


< 0 . 0001 . 


(25) 


Corresponding relations for other cases are obtained from equations by cyclic permutation^. In the cal¬ 

culation we have used the Runge-Kutta method. The accuracy of the calculation is really very good, especially bearing in 
mind the relative simplicity of the equations. The integration was performed with relative precision of 10“®, and the corre¬ 
sponding variant with the precision of 10“^ was exactly the same. One plot contains several thousands of time-steps. Detailed 
information about the initial parameters for all variant s of the calculation is giv en in Table 1. 

The case of a spheroid (a = 6 7 ^ c) is considered in iBisnovatvi-KogarJ (I 2 OO 4 II . Our 3D equations give the same results for 
spheroidal initial conditions. In the case of a three-axis ellipsoid, there are new qualitative effects, because there is an additional 
degree of freedom compared to the case of a spheroid. For the investigation of dynamical behaviour, we start the simulation 
from a spherical body of unit mass (rriin = 1), and zero or small entropy e <C 1. We also specify the parameters characterizing 
different dissipations. In all calculations with relaxation, we use following values: a^ei = 3 and Qei = = ctMi = 15. We 

set non-zero initial velocities. Because of rotation around c-axis, the initial sphere transforms to a spheroid during collapse if 
the initial parameters for axe s a and b are exactly eq ual (a — h, a = b). Note that during the motion spheroids may be not 
only oblate but also prolate llBisnovatvi-Koganll2nn4h . To study the appearance of three-axis figures, the initial parameters 
for axes a and b were slightly different in all variants of the calculations. 

In the first variant (see Fig. 1) there is a large initial angular momentum Min = 0.5. The field of velocities is slightly 
perturbed by increasing b in comparison with a. Initially we observe the collapse and the formation of the pancake. During 
the motion the difference (a — b) increases due to the development of the secular instability, and we obtain the transformation 
of a Maclaurin spheroid into a Jacobi ellipsoid in the dynamics. Because of the relaxation, the oscillation motion damps, 
and the configuration reaches the equilibrium state of the three-axis ellipsoid. The corresponding behaviours of total energy 
Utot, entropy function e and mass m are represented in Figs 2, 3 and 4. We see that the main changes of these parameters 
take place at the stages of the first few oscillations: the entropy function e increases, the mass m and the total energy Utot 
decreases. Than these parameters reach equilibrium values. 


^ Note that expansions given in (42) of|BisnovatiTKogan| contain misprints. The correct expansions at |fc — 1| <g 1 are 


Cfc — 1 + ■ 


1 - fc 


2 2 14 

alfe =- 7(1 - fc); a2fc = --I-—(1 - fc). 

3 5 3 15 


These relations are valid both for A: > 1 and k < 1. 
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Figure 1. The development of an instability at large angular momentum, and the formation of a stationary triaxial figure in case 1 (see 
Table 1 for cases). 



Figure 2. The total energy evolution for case 1. 



Figure 3. The entropy function evolution for case 1. 
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Figure 4. The mass evolution for case 1. 



time 


Figure 5. The relaxation to the spheroid at small angular momentum, case 2. 


In the second variant of the calculations (see Fig. 5) we set a small initial angular momentum Min = 0.1. After the 
collapse, temporary oblate and prolate spheroids appear during the relaxation. There is no secular instability, and the system 
reaches the equilibrium oblate spheroid of Maclaurin. An initial small difference of paramet ers connected with axes a and 
b leads to a small deviation from the results of the similar variant of the 2D simulations JBisnovatvi-Kog^ l2004li . Here 
we have an additional instability: there is an increase in the difference (a — b) at the early stages of motion. However, this 
difference re mains too smal l to be visualized in the figure. T his is the instability characteristic to the system with purely radial 
trajectories JAntonovlll97,?l : rFridman fc Polvachenkolll98R^ . This instability takes place in all variants of the calculations, but 
in most cases it is very small and always disappears during the relaxation and formation of the equilibrium figure. In the cases 
of large angular momentum, this instability does not show up because the difference (a — b) quickly increases as a result of 
secular instability. 

To illustrate the radial instability, we consider an initial sphere without rotation {Mi„ = 0) and without any relaxation 
process. Only the three equations itnii-itTTl with Tnei = 00 have been used, at constant m = 1, M = 0 and e. We have set 
small initial differences between the axes a, b and c, and have obtained an unstable behaviour evidently visible in Fig. 6. 

For numerical investigation of the secular instability, we consider equations of motion without dissipation, and take the 
equilibrium spheroids with small perturbation as the initial configuration. It is easy to obtain the equilibrium parameters 
of the spheroid from the equations of motion with zero accelerations. The initial equilibrium configuration of the Maclaurin 
spheroid was slightly perturbed by increasing b, as in the last column of Table 1. At small angular momentum we obtain 
the conservation of the initial configuration. There are small oscillations around the equilibrium spheroid, during which the 
difference (a — 6) changes sign periodically. At large angular momentum we obtain the development of the secular instability, 
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Figure 6. Radial instability development in the non-rotating body, without dissipation, case 3. 



Figure 7. Development of the secular instability in a rapidly rotating body, without dissipation, case 4. 


and, after several low-amplitude oscillations, the Maclaurin spheroid, is transformed into the Jacobi ellipsoid (see Fig. 7 where 
the variant with large angular momentum is represented). Because of the absence of relaxation, the system does not reach 
an equilibrium conhguration and oscillations continue. If we include relaxation in this variant, we should observe the final 
equilibrium state of the ellipsoid. Our calculations allow us to find approximately numerically the point of the onset of the 
secular instability of a compressible Maclaurin spheroid. It fits well with the analytical results obtained in the next section. 


6 EQUILIBRIUM CONFIGURATIONS AND STABILITY 

Equilibrium of a uniformly rotating figure (spheroid or ellipsoid) is found from the equations 11511 - 1171 with zero time 
derivatives: 

OO 

3m j" du 10 1 £ 25M^ a , . 

2 J (a^ + u)A 3m a (afec)^/® m^ (o^ -I- 6^)^ ’ 

0 



du 

(6^ -I- m)A 


10 1 £ 

3m b (abc)'^^^ 


25M‘^ 
H- — 


b 


( 27 ) 


0 


(a2-p62)2 ’ 
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3m 

0 = -—c 


du 


+ 


10 1 


(c^ + m)A 3m c (abc)^/^ 


(28) 


Introducing the variables 


u , c 
a: = fe = 

a 

and integrals 


/i(fc,fci) = 

73(fc,fci) = 


fei = 


M 
m ’ 


dx 


lo {x + 1)3/2 ^_|_ 2 ;) 

/•OO 7 

/ ax 


/2(fe,fel) = 


da; 


0 (a: + fef)3/2y^(i -(- a:)(fe2 + a;) 


(29) 


(30) 


0 {X + fe2)3/2^(']^ _|_ x){kj + x) 
we can write equations using itTTHi . in the form 


„ 3m r , s 10 £ 

0 = -—a/i(fc,fci)+-^^^ 

„ 3m 10 e 

0 = -—ak^Hk,h)+- 

„ 3m , r r, , \ 10 e 

0 = -—afe/3(fc,fcr)+-^p^. 

Excluding m, we obtain 
Q ^ 1 Ji(fc,fci) ^ ^ ^ 15f (fcfci)^/^ 


25j2 


(l + fc 2 ) 2 ’ 

, ‘25f fci 


(l + fc?)2 ’ 


(31) 

(32) 

(33) 

k^l3{k,ki) ' ^ ' 2e (l + fc2)2’ k^l3{k,ki) ' " ' 2e (1 + fc2)2' 

At g = 6 , fei = 1 equations (E3J are identical, and determine the equilibrium of the Maclaurin spheroid (see iBisnovatvi-KogarJ 
^Oo3)- For Jacobi ellipsoids with a ^ b ^ c, we obtain the following relation between k and fei: 

1 h(k,ki) 


^ ^ kl hjkM) ^ ^ ^ ^2^5 f (fefci)^/® 


Fik,k3) = l-^+ ^ 


= 0 . 


1 D hiKki) k^ l3(k,ki) 

This equation has a trivial solution fci = 1 at all k, corresponding to the Maclairin spheroid. Let us find the bifurcation point 
of the equation 13511 . at which non-trivial solutions appear. While fei = 1 is always a root of equation 1351 . we may write 
F{k, fci) = (fci — l)/(fc, fci). An additional root of the equation F{k, fei) = 0 appears when the root of the equation f{k, ki) = 0 
appears at ki = 1. The root of the zero derivative equation {k, fei) = f{k, fci) + (fei — 1)/^^ (fe, fci) = 0 at fei = 1 coincides 
with the root of the equation f{k, 1) = 0.^ Therefore the value of fc at the bifurcation point is determined by the equation 


dF{k,ki) 

dki 


fei=i = 0 . 


Using and 1351 . this equation is written in the form 

o,.2w,. M , a/2(fc,fcl), a7l(fc,fel), 

2fc 73(fc, fci)|fci=i H--|fei=i-- |fci=i=0 

At fei = 1 we have analytic expressions 


97i(fe,fci), 


=1 = —7o, 


a72(fc,fci) 


afei afei 

so that the equation lEZl is reduced to 


ki=i = —37o, 7o = 


f 


dx 


(1 + x)^Vfe^ + X 


k^ h = Iq, where I 3 = hik, 1 ) and I 2 = hik, 1 ) = h = Ii{k, 1 ). 


Taking the integrals lo, Ii ~ I 2 and I 3 analytically, we have 

fc arccos fe 


73 = ^:^!!- 

Using iunt in iT^ we obtain the equation 

arccosfc fe(13 — lOfe^) 

^ 3 + 8fc2 - 8fe4 ’ 


VI - fe 2 


h=h = - 


fe(l - fc 2 ) 


,2 fe arccos fe 
fe — 


VI - fe 2 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 


We are grateful to A.I.Neishtadt for useful discussion of this point 
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Figure 8. Dependence of the squared eigenfrequency on the axis ratio k in the Maclaurin spheroid without dissipation terms in the 
approximate equations of motion. 


for which the solution k = 0.582724 determines the bifurcation point at the sequence of the Maclaurin spheroids. For the 
uniform spheroid, the position of this point does not depend on the adiabatic index of the matter. 

Above, we have obtained the bifurcation point on the equilibrium curve of the Maclaurin spheroids using only the 
equilibrium relations for the Jacobi ellipsoids. The usual way of investigating stability is to solve linearized equations of 
motion and to find the eigenfrequencies. The linearized equations itTKt-itTTIi without violent relaxation reduce, with account 
of I29II . to the following form: 


5{a — 'h)= 3^ 


-r 


du 


+ 


du 


J (1 + J (1 + 

0 0 


S{a — b). 


The £, Utot, m and M losses are quadratic to perturbations, so these values remain constant in linear approximation. Taking 
5{a — fe) ~ , we come to the characteristic equation 


= 3 ^ 


du 


(1 + u)(fe2 + J (1 + + w)i/' 


OO 

/■ 


du 


(42) 


The plot ud{k) at m = a = 1 is given in Fig. 8. Using the notations of 118811 and lOU in GUi, we write the characteristic 
equation in the form 

u;" = 3^(fe"/3-/o). (43) 

We see that a spheroid loses its stability for the transformation into a three-axial ellipsoid, at the bifurcation point 
k = 0.58272, e = jo? = 0.81267. By derivation it is clear that the position of this point does not depend on the 

polytropic exponent n. Our approximate equations, even without dissipation, describe uniformly rotating ellipsoids (spheroids), 
which are not connected by adiabatic relations, and contain a ’’hidden” non-conservation of the local angular momentum, 
which preserves the uniformity. Therefore, in the presence of this ’’hidden” non-conservation, the loss of stability takes place 
exactly at the bifurcation point. The loss of stability is not connected with relaxation and follows from our equations even in 
the case when all the global integrals r emain constant. In th e exact approach, the instability at this point happens only in 
the direct presence of dissipative terms llChandrasekhaij|l96!jl . The pure adiabatic spheroidal system preserves stability until 
e = 0.952887, where it becomes unstable via a vibrational mode. 

In our approach we cannot exactly investigate the effect of pure bulk viscosity on the stability of the compressible 
spheroids, because the transition is influenced by the hidden non-conservation. Taking account of bulk viscosity and violent 
relaxation does not change the bifurcation point. However, in some cases the presence of relaxation has a stabilizing influence 
on the motion of the system. 

We have considered an equilibrium conhguration of the Maclaurin spheroid close to the bifurcation point, in the interval 
of the secular instability. After perturbation of the density, the spheroid starts to oscillate. Our calculations have shown that 
in the absence of relaxation the development of instability and the transition of the spheroid to a three-axis object both occur. 
In the presence of relaxation the perturbed spheroid remains a spheroid and the system reaches the equilibrium state. This 
behaviour takes place because during relaxation the thermal energy increases, and the eccentricity e decreases and comes into 
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the interval of stability. So the spheroid initially perturbed in the interval of instability, not far from the bifurcation point, 
becomes stable due to the p rocesses of relaxation. 

According hypothesis of lOstriker fc PeeblesI Jl97,^l , the stability of an isolated axially symmetric system is determined by 
the ratio Urot/ \Ua\- They determined from numerical experiments the critical value for various configuration as 0.14 ± 0.03. 

found that compressible spheroids become secularly unstable to triaxial deformations at the bifurcation point. 


Eh 




where Urot/\Ug\ — 0.1375, independent of n. Our formula gives exactly the same result, which is also confirmed by our 
numerical simulations. As noted above, this result remains valid also in the presence of dissipation. 


7 CONCLUSIONS 

We have investigated the dynamics of a three-axis dark matter ellipsoid. The equations of motion for the axes of a uniform 
compressible ellipsoid have been obtained by variation of the Lagrange function, in which violent relaxation and losses of 
matter, energy and angular momentum have been included phenomenologically. 

The system was solved numerically, until the formation of stationary rotating figures in the presence of relaxation. For 
lower angular momentum M we have the formation of a compressed spheroid, while at larger M we follow the development 
of a three-axial instability and the formation of a three-axial ellipsoid. The instability in this approximation happens at the 
bifurcation point of the sequence of Maclaurin spheroids, where the Jacobi ellipsoidal system starts. 

The bifurcation point coinciding with the point of loss of stability is found analytically in the form of a simple formula, 
by static and dynamic approaches. Numerical and analytical considerations give identical results. 

The development of instability, connected with radial orbits, is obtained for slowly rotating collapsing bodies. 
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